epi draft 



Non-equilibrium length in granular fluids: 
|from experiment to fluctuating hydrodynamics 



o : 

(N ■ 

Q_G. Gradenigo, a. Sarracino, D. Villamaina and A. 
CNR-ISC and Dipartimento di Fisica, Universitd Sapienza 

o ■ 
(N : 



PUGLISI 

- p.le A. Mora 2, 00185, Roma, Italy 



PACS 45 . 70 . -n - Granular systems 
PACS 05 . 20 . Dd - Kinetic theory 

PACS 05. 40. -a - Fluctuation phenomena, random processes, noise, and Brownian motion 

Abstract. - Velocity correlations in a quasi-2D driven granular fluid are studied in experiments 
and numerical simulations. The transverse velocity structure factor reveals two well deflned energy 
scales, associated with the external "bath temperature" Tt, and with the internal granular one, 
Tg < Tt, relevant at large and small wavelengths respectively. Experimental and numerical data 
are discussed within a fluctuating hydrodynamics model, which allows one to define and measure 
a non-equilibrium coherence length, growing with density, that characterizes order in the velocity 
field. 
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The study of collective phenomena in complex systems 
is a central issue in statistical mechanics. In particular, 
the degree of order emerging in such systems is customar- 
ily investigated by studying the behaviour of correlation 
i\inctions [1]. However, in many systems which display 
interesting collective phenomena interactions are non- 
tonservative and a stationary state is achieved through a 
continuous and homogeneous energy injection acting on 
every single particle: this is the case, for instance, of 
homogeneously vibrated granular materials as well as of 
the so-called "active fluids" including swarms of bacteria, 
birds, insects, fishes, pedestrians and even artificial mo- 
tors [2-4]. In all those systems, the study of correlations 
requires tools beyond the standard equilibrium ones. 
. The emergence of order out of equilibrium is studied 
here in a granular system [5,6], where inelastic collisions 
induce cross correlations between particles velocities. The 
behaviour of such correlations is customarily described 
with fluctuating hydrodynamics (FHD), which relies on 
a separation of time and length scales in the system, to- 
gether with a suitable model for noise. Assessing the limits 
of validity of granular hydrodynamics is non-trivial [7, 8] , 
but in several situations it has proven to be useful in ex- 
plaining the behaviour of large scale correlations in gran- 
ular fluids [9, 10]. In a more general context, the existence 
of non-trivial correlations emerging from non-equilibrium 
assumptions was already noticed in [11]. Recently, FHD 
also provided a prediction for the slowing-down of the dy- 
namics, with increasing packing fraction, in granular flu- 



ids [12, 13], a phenomenon which has a well known coun- 
terpart in elastic systems [14]. The above mentioned and 
similar results [15], which clearly show the existence of a 
time-scale growing with the density, suggest to look also 
for a growing length-scale. 

In this paper we show the existence of a finite correlation 
length ^ in the velocity field, peculiar to a non-conservative 
fluid. Velocities correlations are investigated studying the 
velocity structure factor measured from experiments and 
numerical simulations. Data are discussed within a FHD 
model, allowing us to extract a length scale ^ and conflrm- 
ing that its associated order is a peculiarity of the out-of- 
equilibrium regime. The behaviour of ^ is investigated at 
different packing fractions, showing a remarkable increase 
of ordering in the velocity fleld with growing density. The 
large scale fluctuations of the velocity fleld suggest that 
a reference theoretical model where the thermostat has a 
finite viscous drag [16] is the most appropriate to modelize 
the experiments. 

The experimental setup, sketched in Fig. 1, consists of 
an electrodynamic shaker (LDS V450) which vibrates ver- 
tically (£ axis) a horizontal rigid aluminium plate cov- 
ered by a partial granular monolayer of N "steel 316" 
spheres, with packing fraction (p = -kuo^ /A < 0.5, where 
n = N /Area is the number density and the diameter a of 
spheres is 4 mm. The plate is circular and has slowly rising 
boundaries in order to prevent highly dissipative head-on 
collisions between spheres and walls. The radius of the 
central flat region is i? = 100 mm. The surface of the plate 
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is rough [17-19], transferring part of the vertical momen- 
tum to the horizontal (xy plane) motion of the spheres. 
The plate is roughened by stainless-steel sand-blasting 
with a granulometry in the range [100— 500]/xm. Collisions 
among the spheres are also responsible for z ^ xy energy 
conversion. From a fast camera (Mikrotron EoSens CL, 
up to 506fps for 1280x1024 pixels), placed on the top of 
the plate, we observe the xy motion of the spheres. Even 
if the container has no lid (as it was in the experiment 
with the rough plate in [17]), we verified by visual inspec- 
tion that - at the chosen shaking parameters (see below) - 
particles never overcome each other and constitute a quasi- 
two-dimensional fluid. The results discussed below are ob- 




Fig. 1: (color). Schematic of the experimental setup. 



tained for sinusoidal shaking, i.e. the plate follows the law 
z{t) = Asin(wt), with Auj"^ / g = 12 and w/(27r) = 150 Hz, 
where g is the gravity acceleration, yielding a maximum 
velocity of the vessel Vm — 125 mm/s. A particle track- 
ing software identifies particles' centers of mass, up to an 
accuracy of 1/5 of a pixel, and their displacement along 
close successive frames (Ai = 1/57 sec), reconstructing 
particles' velocities: the analysis of the mean squared dis- 
placement As^(r) = {\x{t -\- t) — x{t)\'^) guarantees that 
the reconstruction is operated in the ballistic regime, i.e. 
As2(t) ~ for T < 2Ai. 

The measure of the static structure factor Sp{k) = 
1/y (I exp(— ik • r.i)p), where r.; denotes the posi- 
tion of each particles and V is the volume, reveals that 
in the range of packing fraction that we studied, i.e. (p € 
[0.1, 0.42], the structure of the granular fluid presents small 
differences from its elastic counterpart (inset of Fig. 2), 
probably due to weak clustering or inhomogeneities. The 
small enhancement of the structure may also be under- 
stood by realizing that - in our setup - a slight amount of 
vertical motion allows particle to arrange at separations 
even smaller than their diameter. 

More peculiar is the behaviour of the velocity structure 
factor. In particular we study here its transverse part: 

n^S^k) = V-^v^{k)v^{-k)), (1) 

with v±{k) — J2i (^^i ■ ^'-l) cxp(— zk-r^), where Vj denotes 
the velocity of the i-th particle and k±_ is the unitary vector 
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Fig. 2: (color online). Sj_{k) measured in units of in 
experiments (where Vm is the maximum velocity of the vi- 
brating plate), for different values of (f). Data represent av- 
erages over 3300 pairs of frames. The dotted lines on the right 
show the corresponding Tg. As a reference for the eye, we 
also indicate the average (among the four experiments) of Tt, 
which is defined in the text. Inset: Sp{k) for (j) = 0.14 (a) 
and (f) — 0.42 (•) with the analytical predictions for elastic 
disks [20]. In all figures the lowest wave- vector is ka > 2kminO 
{2k,nin<y ~ 4-Ka/R = 0.5), because at smaller ones finite size- 
effect are dominant. The values of S±{k) reported are an av- 
erage on moments with modulus k £ [n2n/R, {n+ l)2n/R] for 
integer n > 2. 



such that kj^ ■ k ~ 0. The curves of nSj_{k) measured in 
the range of packing fractions (j) £ [0.1; 0.42], i.e. in a 
moderate dense fluid regime, are shown in the main frame 
of Fig. 2. 

By definition, S±{k) yields a measure of the energy lo- 
calized on average on each mode k. In the granular system 
studied here, we find that a different amount of energy is 
concentrated on each mode, namely at stationarity S±{k) 
is a non trivial function of k. More specifically, we observe 
that two relevant energy scales are well defined in the sys- 
tem: one at large fc, which corresponds to the granular 
temperature Tg = i/{'2N)^-{vf) and strongly depends 
on the packing fraction (j); the other, at small fc, is weakly 
dependent on 0. Consider that in the limit fc — > 0, we 
have that, by definition, nS±{k) — > N{V^^), with Vcm the 
velocity of the center of mass of the system. The quan- 
tity N{V^^) is not affected by momentum-conserving col- 
lisions: it is, therefore, related only to the interaction with 
the vibrating plate. For such a reason, we use this value 
as a measurement of the so-called "bath temperature" Tt,. 
In the theoretical model introduced below this energy cor- 
responds to the temperature of the thermostat coupled to 
the granular fluid. From Fig. 2 it is evident that the gap 
between the two energy scales defined above significantly 
grows as the packing fraction is increased. We will show 
that the existence of such a gap is deeply connected with 
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the existence of a finite correlation length in the veloc- 
ity field. The oscillations of S±{k) at high values of k, are 
likely to be a signature of an unexpected coupling with the 
density modes, already observed in previous works [21]. 

With the aim of giving a consistent interpretation of our 
experimental results, we now propose an effective model 
for the microscopic dynamics of particles (restricted to 
the plane xy) and deduce from it a theoretical formula 
- through fluctuating hydrodynamics [21,22] - for S±{k). 
Our proposal is well known in the literature [16]: it mim- 
ics the action of the rough vibrating plate by means of an 
interaction with an effective viscous "bath" . The model is 
described by the following equation governing the dynam- 
ics of the i-th particle (i S [1, A^]) 

Mt)^~lbMt)+C^{t)+F,, (2) 

where represents the effect of particle-particle inelastic 
hard core collisions with restitution coefficient a, while 
7(^^ is the typical interaction time with the plate and 
Ci(t) a white noise, with zero mean and (Ci,a(i)Ci,^(^')) = 
2Tbjb5ci3S{t - t')5,j (a,^ = {x,y}). This kind of ther- 
mostat is well defined also in the elastic limit a — >■ 1, 
where the fluid equilibrates to the bath temperature: 

For the above model one can apply the usual kinetic the- 
ory to study hydrodynamic fields [21]. In linearized hydro- 
dynamics, which applies to the moderately dense regime 
here studied, the transverse modes v_\_{k, t) of the velocity 
field are decoupled from all the other modes. In particu- 
lar, at small k they are well described - on average - by 
the equation v±{k,t) ~ —(76 + h'k'^)v±{k,t), with v the 
kinematic viscosity [1]. This equation contains two main 
sources of relaxation: the external friction parametrized 
by 7b, and vorticity diffusion parametrized by i^. To an- 
alyze structure factors one needs to plug in appropriate 
noises. The two relaxation mechanisms described above 
are associated with two independent sources of random- 
ness: external and internal noise, respectively. A rigorous 
evaluation of the way those noises enter the description for 
v±{k,t) can be found in the literature for different mod- 
els [22]. Here we present a more phenomenological ap- 
proach, assessed by its ability to reproduce experimental 
and numerical results. The idea, partly inspired to previ- 
ous granular FHD theories [21], consists in keeping valid - 
as a first approximation - the 2nd kind Fluctuation Dissi- 
pation Relation (FDR) for the two relaxation-noise pairs, 
separately. The result can be cast in a Langevin equation 
for the shear mode, with one single noise with appropriate 
variance: 

v±{k,t) = -{■jb + iyk^)v±{k,t)+r]{k,t), 
{■n{k,t)T^{-kX)) = 2{-ibn + vk^Tg)5{t-t'), (3) 

with {'q{k, t)) = 0. The amplitudes of the external and in- 
ternal noises in the above equation are reasonably assumed 
to be proportional to Tt, and Tg, respectively. The p.d.f. 



of vj_{k) measured in experiments are Gaussian, in agree- 
ment with Eq. (3), whereas those of real space velocities Vi 
show deviations from Gaussianity due to the relevance of 
collisions in Eq. (2). From Eq. (3) wc immediately obtain 

nSAk)^N-\\vAkr)^:^^I^±^. (4) 

7h + vk'' 

At equilibrium, namely when collisions are elastic and 
Tg = Tf,, cquipartition between modes is perfectly fulfilled 
and the structure factor becomes flat, i.e. Sj_{k) ~ Tt. 
Differently, in the granular case, where Tg ^ Tb, cquipar- 
tition breaks down and from Eq. (4) we have that S± — >■ Tf, 
for small k and S± — > Tg for large k. The Fourier inversion 
of Eq. (4) in two dimensions yields: 

nG^(r) = r,<5(^)(r) + (T, - Tg)^^^, (5) 

where Ko{x) is the 2nd kind modified Bessel function that, 
for large distances, decays exponentially 

and we have introduced ^ = v/^b- Such quantity turns 
out to be a non- equilibrium correlation length, measurable 
in the system only when non-conservative interactions are 
present. Indeed, at equilibrium, the second term in the 
r.h.s. of Eq. (5) vanishes and no coherence can be identi- 
fied in the velocity field fluctuations. 

The saturation of nS±{k) to a flnitc value at small k pro- 
vided from this theory is in agreement with the behaviour 
of experimental data observed in Fig. 2, suggesting that 
the model with finite friction is appropriate to describe 
the system. 

The correlation length <^ is defined through 7t and ly, 
which can be obtained by studying the dynamical cor- 
relation Cj_(A:, ^ (w±(A:,t)?;j_(-fc,0))/(|wj^(fc,0)|2). In- 
deed, from Eq. (3), we immediately obtain C±{k,t) = 
exp{— (7(, -I- vk^)t}. Such behaviour is in good agree- 
ment with experimental data, as reported in Fig. 3, where 
C±{k,t) is plotted for several values of (Fig 3(a)) and 
different wave- vectors k (Fig. 3(b)), together with the best 
fit curves. At each packing fraction, the length scale 
C = \/Wjb can be obtained by comparing the decay of 
shear modes at different k (see Fig. 3(d)). A parabolic fit 
T~^{k) = jb + i^k'^ shows that our data are broadly con- 
sistent with the theory, but it is also important to notice 
that only few data points are available. In Fig. 3(c) we 
see how the correlation length ^, rescaled with the pre- 
diction - in the diluted approximation - for the mean free 
path Ao = l/(20F(?2nCT), where 52 = (l-7(/./16)/(l-(A)2, 
significantly grows with the packing fraction and signals 
an increasing degree of order in the granular fluid. No- 
tice that the main contribution to the growth presented 
in Fig. 3(c) is due to the decrease of Ao, while ^ grows 
only weakly, as already visible in Fig. 2, where S±{k) has 
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Fig. 3: (color online). Data from experiments, (a) Temporal 
decay of C±{k,t) at fixed k = 2kmin and (b) for different k 
at fixed (j> ~ 0.42. The unit of times is 2tt/ijj. (c) Correlation 
length 5, rescaled by Ao, plotted as function of (j>. (d) Parabolic 
fit of r-i(fc) =7i, + !yA;^ 



a "sigmoidal" shape with an inflection point almost inde- 
pendent of 4>. Such an observation is consistent with the 
experimental measurements discussed in [17]. The growth 
of ^/Ao denotes an increasing number of correlated parti- 
cles. 

The finite extent of velocity correlations in a granular 
fluid was already pointed out in [17] . Here we discover that 
it is strictly related to a saturation of S±{k) as fc — >■ 0, 
corresponding to a well defined energy of the "external 
thermostat" . Within the homogeneous cooling regime a 
length-scale ^ was discussed in [10], but in that context it 
separated the stable, k > from the unstable k < 
modulations of the velocity field. 

The above discussion shows that a fluctuation on the 
mode k decays with characteristic time T(fc) = l/(7ti -I- 
i^fc^). This timescale has to be much larger than the 
microscopic time scale given by the mean collision time 
Tc ~ ^o/\/T^, for the FHD with white noise to be effec- 
tive. Comparing T(fc) and Tc, we see that in our experi- 
ments such constraint is fulfilled for a narrow interval of 
small k. This range of validity is constrained by parame- 
ters such as the kinematic viscosity v and the drag of 
the thermal bath, which are not under control in experi- 
ments. 

The predictions of the linearized FHD described above 
arc tested by simulating an interacting system of inelastic 
hard disk modclizcd by Eq. (2), where 7^ enters as a fixed 
parameter. In particular, 7^^ ^ Tb is the timescale of the 
thermal bath, which can be tuned in order to have smaller 
values of Tg. Indeed, by lowering 7^, also v is lowered. 



r(/c) = l/(7b -|- I'k^), extending the range of validity of 
FHD, i.e. the interval of k values where Eqs. (3) and (4) 
are expected to hold. We have performed 2D noisy event- 
driven molecular dynamics simulations of inelastic hard 
disks, at different packing fractions cj) G [0.1,0.5], with 
a e [0.6,1]. Details and results of the simulations are 
reported in a separate publication [24]. We consider N = 
10000 disks in a box of area with periodic boundary 
conditions and the packing fraction is varied by changing 
the side length of the box. The particles have diameter 
a = 0.01 and unit mass and the thermal bath of Eq. (2) 
has parameters = 1 and 7;, = 1, while Tc is always in 
the range [0.005,0.04]. In Fig. 4 we report our results for 
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because, in first approximation, u ^ Tg^^ [23]. Therefore, 
by tuning 7^, we can raise the shear mode decay timescale 



Fig. 4: (color online). S±{k) measured in event driven molecu- 
lar dynamics simulations, averaged over 400 frames. The solid 
lines are the best fit curves via Eq. (4), for = 0.42 and 
(j) — 0.14. Inset: (,/Xo as a function of 0. 

S±{k). A striking analogy with the experimental results 
of Fig. 2 can be appreciated. Again, two energy scales are 
clearly observed: at small k all curves saturate to the value 
T;, = 1, in agreement with Eq. (4). Moreover, as expected, 
for a = 1 energy equipartition holds and S±{k) = Ti,, 
whereas for a ^ 1 the structure factor at large k oscillates 
around the granular temperature Tg. Notice that here, 
by using Eq. (4), we obtained a good fit of data, in a 
wide range of k, thus validating the FHD description (solid 
lines in Fig. 4). The kinematic viscosity v obtained from 
these fits, yields the correlation length ^. A remarkable 
increase of ^/Aq as function of is found (see the inset of 
Fig. (4)). Notice also that the gap between the granular 
temperatures at different packing fractions and the bath 
temperature is strongly increased, compared to the one 
observed in experiments. Such a large separation between 
the two energy scales is consistent with an extended realm 
of validity of the FHD. 

In conclusion we have measured the transverse velocity 
structure factor, and its time decay, for a 2D driven gran- 
ular fluid. We have observed a striking agreement with 
simulations of a microscopic model with stochastic driving 
and viscous friction, and with a fluctuating hydrodynam- 
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ics theory with two noise sources at different temperatures 
Tb > Tg. A central issue of our study has been to show 
the increasing extent of correlations in the velocity field, 
with growing density. Such a phenomenon emerges in the 
system only when time-reversal symmetry is broken, i.e. 
in the presence of an energy flux from the external driving 
to the system, due to the inequality Tg < Tb- 
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